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Understanding the thermo-mechanical response of the Space Shuttle External Tank spray- 
on foam insulation (SOFI) material is critical to NASA’s Return to Flight effort. This closed- 
cell rigid polymeric foam is used to insulate the metallic Space Shuttle External Tank, which 
is at cryogenic temperatures immediately prior to and during lift off. The shedding of the 
SOFI during ascent led to the loss of the Columbia, and eliminating/minimizing foam loss 
from the tank has become a priority for NASA as it seeks to resiune scheduled space shuttle 
missions. Determining the nature of the SOFI material behavior in response to both thermal 
and mechanical loading plays an important role as any structural modeling of the shedding 
phenomenon is predicated on knowledge of the constitutive behavior of the foam. 

In this paper, the SOFI material has been analyzed using the High-Fidelity Generalized 
Method of Cells (HFGMC) micromechanics model, which has recently been extended to 
admit a triply-periodic 3-D repeating unit cell (RUC). Additional theoretical extensions that 
were made in order to enable modeling of the closed-cell foam material include the ability to 
represent internal boundaries within the RUC (to simulated internal pores) and the ability 
to impose an internal pressure within the simulated pores. This latter extension is crucial as 
two sources contribute to significant internal pressure changes within the SOFI pores. First, 
gas trapped in the pores during the spray process will expand or contract due to temperature 
changes. Second, the pore pressure will increase due to outgassing of water and other species 
present in the foam skeleton polymer material. With HFGMC’s new pore pressure modeling 
capabilities, a nonlinear pressure change within the simulated pore can be imposed that 
accounts for both of these sources, in addition to standard-thermal and mechanical loading; 

The triply-periodic HFGMC micromechanics model described above was implemented 
within NASA GRC’s MAC/GMC software package, giving the model access to a range of 
nonlinear constitutive models for the polymeric foam skeleton material. A repeating unit cell 
architecture was constructed that, while relatively simple, still accounts for the geometric 
anisotropy of the porous foam microstructure and its thin walls and thicker edges. With the 
lack of reliable polymeric foam skeleton material properties, many simulations were executed 
aimed at backing out these material properties. Then, using these properties, predictions of 
the thermo-mechanical behavior of the foam, including calculated internal applied pressure 
profiles, were performed and compared with appropriate experimental data. 
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Multi-Scale Approach to Modeling ET Foam Insulation 


Global Tank FEA 


Standard Elements 


Line 
Elements 


Integratio n Point Level HF GMC Analysis of Foam Material 


large Defects Modeled 
Explicitly in FEA 


Note : Each element can be assigned 
its own geometry, material 
properties, thermal history, and 
internal pressure history 


Local Foam FEA 
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NASA/OAI ImMAC Software Suite for 
Multi-Scale Analysis of Advanced Materials & Structures 
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I Hyper-MAC called from within HyperSizer 
] FEAMAC called from within ABAQUS 
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HFGMC Micromechanics Model: Geometry and Approach 



HFGMC models a periodic 

material 

Repeating Unit Cell (RUC) is 
identified 

RUC is discretized into an arbitrary 
number of subcells (apy) 

Subcells may contain any material, 
or be empty (pore) 

New capabilities included: 

Internal Pore Pressure 
Elastic Material Stiffening 
Newtonian-Viscous Material Creep 


HFGMC Micromechanics Model: Basic Equations 


• Empty subcells subjected to internal pressure: F ^ r) 

* Constitutive equation within each subcell including inelastic and thermal strains: 

= c («fr) / z {<&) _ z n<#r) \ _ j**r) A T 


• Second-order displacement field assumed in each subcell: 

u<^> =*x+ + 




• Unknown terms in each subcell, determined via imposition of equilibrium, 
continuity, and periodicity equations in an average (integral) sense 


• Results in system of linear algebraic equations for unknown terms: KU = f + g 

K - Contains geometry and thermo-mechanical property information 
U - Contains unknown terms ( ) 

f - Contains applied strain and temperature information 

g - Contains integrals of inelastic strains and effects of system of pressures, P {afir) 


HFGMC Micromechanics Model: Basic Equations 

Solving the system of equations establishes the localization relation: 

z ( ^=A {afir) £ + A th{afir) AT + A ,{afir) + A P{aflr) 

Mechanical Thermal Inelastic New 

concentration concentration concentration Internal pressure 

tensor vector vector concentration vector 


• The global average stress in the material-filled subcells are related to the average 
subcell stresses by: 

*=ik £ 




• Above equations allow establishment of constitutive equation for porous material: 

a = C*e-(rA7' + a , )-a' > 


Where, 
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Modeling Smooth Tensile Test 
Specimen Results 

LN2 Cycle following RT Load 
BX265, Manual Spray, 1-1 Orientation 
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Implementing Stiffening and Creep 

• Models need to be relatively simple, but capture primary effects 

• Elastic Stiffening Constitutive Model (Chen and Saleeb, 1981): 

- Bulk Modulus: = /<"„ + K x e v + K 2 e*. 

volumetric strain: £„ = £,,+ e 22 + e, 3 

- Shear Modulus: g(p,^7J)=C 0 +^/?+^ > /^ 

P=j(oii+<* 21+*3») 

- For simplicity reduce to two parameters, set: K 2 = A, = 0 

• Newtonian-Viscous Creep Model for Polymers (Ashby and Jones, 1980): 

£ creep = ^<Te~ G/ * r Q = Creep activation energy 

A = Creep coefficient 


Creep and Stiffening Model for 
Skeleton Material 



Notes 

* Only tensile causes 
stiffening 

•Creep strain does not 
contribute to stiffening 

•Creep highly 
dependent on 
time/loading rate 

•Temp-dependence of 
creep automatic 
through Arrhenius eq, 

•Additional temp- 
dependence of G 0 and 
Kq not shown in chart 


5 


Model Variables 

• Geometry 

- Aspect Ratio 

- Wall Thickness 

- Edge Thickness 

- Unit Cell Discretization 

y Treat as fixed 

• Material Properties 

- Elastic: E 0 (T), v 0 = 0.4 

- Elastic Stiffening: K t , A 2 

- Creep: A, Q 

- Thermal Expansion: a(T) 

"\ 

y Attempt to determine to 
achieve good agreement 






Room Temp Loading and 15 s Creep 


Note 



Cool Down to -320 F 

• Interaction of many effects : 

- E 0 temp-dependence 

- Pore pressure change 

- CTE and CTE temp-dependence 

- Creep: temp-dependence due to Arrhenius term 

• Temp, change rate (significantly affects creep) 


* temp-deoendence : 
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• Utilized E vs. T data 

• -320 F vs. 68 F — * E increase by 46% 

• Assume linear increase in skeleton 
material E 0 from room temp to cryo 

E = £ 0 [l -(0.0021 15 °C~ i )A7’] 

= £ 0 [l -(0.001175 T _, )Ar] 


Cool Down to -320 F 


• Pore pressure change : 

- Assume no phase transfer of volatiles from room temp to cryo 

- Assume ideal gas contraction of air in pores 

ap - PwJi KT (Q-QQ1 \979 g/cc)(83\4J/ mole -K) (lOOcc) 3 _ 
* M air 28.8 g I mole m' 

= ( 345.8 Pa/K)AT = ( 0.0279 psij‘F)AT 


• CTE and CTE Temp-Dependence: 



-«» -WO -}V0 -300 -250 -100 -150 -100 -50 0 S« 100 ISO 


of BX 205 (3*) in che 1 1 dwwnoo micr V kwh and 
A 0.9 pw load 


• Used data from Southern 
Research to back out CTE 
(1 atm data) 
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Prediction of Southern Research 
Vacuum Data 


15 



200 


Exp (vac) 
Exp (vac) 
Exp (vac) 
Model (vac) 


Test Sequence 

• Cool to cryo at latm 

• Pull vacuum 

• Heat to 125 F 

• Cool to 75 

• Release vacuum 

Model Sequence 

oii»t Time (ft) Tcmp(F) Pext(p>i) Putt (pei) 
I 0 75 14.7 0.00 
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9200 
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125 
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- 13.51 

1.19 

14.70 

23.1* 

14.70 

000 


Temperature (F) 


Cool Down to -320 F 


LN2 Cycle following RT Load 
BX265, Manuel Spray, 1-1 Orientation 



* Last item needed is cool down rate (temperature vs. time) 

* First tried linear temp. vs. time 


Cool Down to -320 F 

(assuming linear time vs. temp.) 



Transient Heat Transfer in Cylinder 



Analytical Solution (Bowman, 1958) 


T{r,t) = T b -2(T„- «/* ) 

«=l a n J ]\ a n ) 

J 0 , J. Bessel functions of the 1 st kind 


Evaluate using 
Mathematica 


nth positive zero of J 0 



Transient Heat Transfer in Cylinder 


• Also performed ABAQUS FEA heat transfer analyses: 



Transient Heat Transfer in Cylinder 


ABAQUS Temp Profiles 
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Transient Heat Transfer in Cylinder 


Average Temp vs. Time 



Cool Down to -320 F 



Cool Down to -320 F Inconsistent 
with Southern Research Data 



Re-calibrate CTE 



Full Smooth Tensile Test Simulation 



GRC Thermal Expansion Data (Feb.2006) 

Thermal Expansion of BX-265, Block 3 



Elevated Temperature Pore 
Pressure Estimates 



GRC Thermal Expansion Test Simulation 



Conclusion 

• NASA/OAI ImMAC suite enhanced to enable micromechanics analysis 
of foams with: 

- Internal pore pressure 

- Elastic material stiffening 

- Newtonian-viscous material creep 

• Model attempts to captures 1 st -order effects of many interacting 
mechanisms 

• Model can do a reasonable job of simulating foam material thermo- 
mechanical behavior 

- Quite a few parameters must be backed out of foam test data 

- Micromechanics more powerful when constituent-level data are available 


